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A random matrix approach to glassy physics is introduced. It leads to a class of models which exhibit 
both, glassy low-temperature phases, and double- and single-well configurations in their potential energy. 
The distribution of parameters characterizing the local potential energy configurations can be computed, 
and differ from those assumed in the standard tunneling model and its variants. Still, low-temperature 
anomalies characteristic of amorphous systems are reproduced, and we are able to distinguish properties 
which can be expected to be universal from those which cannot. 
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Glassy materials exhibit, at low temperatures, a num- 
ber of properties which are considered anomalous in com- 
parison to those of their crystalline counterparts . Ex- 
amples are the roughly linear temperature dependence 
of the specific heat CiT) ||] and the approximately 
quadratic temperature variation of the thermal conduc- 
tivity k{T) at T < 1 K, which are in contrast to the 
behaviour of these quantities in crystals. Moreover, be- 
tween approximately 1 and 20 K, one observes a crossover 
to a behaviour of the specific heat and a plateau in 
the thermal conductivity. The anomalies below 1 K ap- 
pear to be universal in the sense that they are shared by 
a large variety of amorphous systems, whereas between 
approximately 1 and 20 K a stronger dependence on prop- 
erties of the specific material appears. At still somewhat 
higher temperatures, a certain degree of universality is 
again observed, at least in thermal conductivity data 

It is the existence of a broad range of localized low 
energy excitations which is generally held responsible for 
the aforesaid anomalies. The dominant mechanism giv- 
ing rise to excitations below 1 K is believed to be tun- 
neling in double- well potentials (DWP). This is the con- 
tent of the phcnomenological standard tunneling model 
(STM) [§,1). A different set of excitations in amorphous 
systems can exist as localized vibrations in soft anhar- 
monic single-well configurations of the potential energy. 
They have been postulated within the, likewise phe- 
nomenological soft -potential model (SPM) to de- 

scribe the physics also above 1 K. Either model requires 
specific assumptions concerning the distributions of its 
parameters to describe glassy physics at low T. 

The existence of localized soft vibrations ||], and 
DWPs H] responsible for two-level tunneling systems 
(TLS) has been convincingly established in molecular 
dynamics studies of glassy systems. In the statis- 
tics of the parameters characterizing the local potential 
energy configurations was also investigated (within the 
confines of a generalized SPM). Still, as of now it is per- 
haps safe to characterize both, STM and SPM as provid- 
ing phcnomenological descriptions, based on assumptions 
which, while plausible in many respects, are still lacking 



support based on more microscopic approaches, such as 
that of for KBr:KCN mixed crystals. At the same 
time they are widely used to explain and fit experimen- 
tal results. In particular, it has been argued that 
the STM could not really account for the considerable 
degree of universality observed in glassy physics, and nu- 
merical investigations on specific substances can at 
best only partially elucidate this point. It was suggested 
]lOt that universality in glassy systems would emerge only 
due to sufficiently long-range interactions between TLSs, 
an idea which has been taken up in [ p^ . 

Recently, one of us has proposed a random-matrix 
model of a glass, which exhibits both, an amorphous 
low-temperature phase, and a collection of double- (and 
single-)well configurations in its zero-temperature po- 
tential energy surface The main motivation for 
that study was to formulate a bona-fide model of an 
amorphous system for which the statistics of local poten- 
tial energy configurations becomes part of the world of 
the computable. The result, taken as input for comput- 
ing the influence of tunnelling-excitations on the low- 
temperature properties, was used to demonstrate that 
the model shares the low-temperature anomalies charac- 
teristic of glassy systems. Moreover, it was shown that 
one could also exhibit relations between low-T and high- 
T phenomena, e.g. between the low-temperature specific 
heat and the value of the glass-transition temperature it- 
self, and relate this, in turn, with features of the system 
at the microscopic level of description. 

The purpose of the present paper is, inter alia, to in- 
vestigate the model beyond the two-level approximation 
in the quantum-mechanical treatment of its local exci- 
tations. Thereby we find that it reproduces — without 
further assumptions — also the characteristic bump as it 
is observed in C{T)/T^ plots of the specific heat. We can 
trace this back to microscopic features of our model, and 
we are able to address the universality question posed in 
]lO| : Our answer, in short, is that the ensemble of local 
potential energies is itself a largely collective affair, and 
some degree of universality may therefore be expected to 
arise without interactions at the level of quantized local 
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excitations. We shall elaborate on these points below, 
and also try to situate our model in glassy physics, be- 
yond the specific assumptions embodied in it. 

Let us restate the main ingredients of our model. It is 
based on the following ansatz for the potential energy of 
A'' degrees of freedom ("particles") forming a glass-like 
system, 

Upot{v) = X! -^y^^^J + r X! ' 

in which Vi may be interpreted as the deviation of the 
i-th particle from some preassigned position. The amor- 
phous aspect is modelled by taking the first, harmonic 
contribution to J7pot to be random, so that the reference 
positions would generally turn out to be unstable at the 
harmonic level of description. A set of anharmonic on- 
site potentials G{vi) is therefore added to stabilize the 
system as a whole. We choose the harmonic part such 
that the model can be analysed exactly within mean-field 
theory, and the replica method is used to deal with the 
disorder. Namely, we take the Jij to be Gaussians with 
mean Jq/N and variance J'^/N, specializing to J = 1 to 
fix the energy-scale. For the on-site potential we choose 

G{v) = lv' + ^v' ■ (2) 

That is, G also creates a harmonic restoring force, and by 
varying 7 we can tune the number of modes in the sys- 
tem which are unstable in the harmonic approximation. 
Other forms of G{v) may be contemplated; our ability to 
solve the model does not depend on the particular shape 
of G. The only requirement is that it increases faster 
than for large \v\ for the system to be stable, and our 
choice is the simplest respecting a Z2 symmetry. 

The choice of random couplings in (1) puts our model 
outside the class of glass-models in the narrow sense. In 
view of recent ideas concerning the fundamental similar- 
ity between quenched disorder and self-induced disorder 
as it is observed in glassy systems proper |l8j , it may still 
be argued that our choice should capture essential aspects 
of glassy physics. Models of this type with different cou- 
plings and different choices for G{v) have been studied 
in the context of analogue neuron systems ||l^,|l^ . 

In a mean-field setting, the system is described by an 
ensemble of effective single-site problems, characterized 
by potentials which, in replica theory, attain random pa- 
rameters. Solving the model means to compute the dis- 
tribution of parameters characterizing the local potential 
energy configurations self-consistently, and this is pre- 
cisely what we were afer: it can directly be translated into 
the distribution of the parameters characterizing DWPs. 
Taking this distribution as input of a tunneling model, 
we were able to compute, for instance, the contribution of 
the tunneling states to the specific heat at low T . In 
what follows, we shall go beyond the two-level tunneling 
approximation, and begin to explore the consequences. 



To analyze the potential energy surface, one computes 
the (configurational) free energy of the system 

fNiP) ^ -iPN)-Hii fUdv.eM-PUpotiv)] (3) 

and takes its T = limit to eliminate entropic contribu- 
tions, and to select one of the system's (many) ground 
state configurations, using replica theory to average over 
the disorder so as to get typical results. Standard argu- 
ments give /(/?) = hm„^o/n(/3) for the quenched 
free energy, with 

n/n(/?) = ^JoEp' + i/5E«'6 

a a,b 

- (^'^ / n '^"^ «=^p [ - /^f^cft ({«"})] • (4) 

Here 

C/cff = -JoY.PaV'' - ^ ^ qabV'^V^ + - ^ G{v'') (5) 
a a.h ^ a 

is a replicated single-site potential and the order param- 
eters Pa = N-^ Y.I i'^i) and qab = N'^ (K'^t) are 
determined from the fixed point equations 

Pa = (O , a - l,...,n (6) 
Qab = , a,6= l,...,n , (7) 

in which (...) denotes a Gibbs average corresponding to 
the replica potential (5), and the limit n — )■ is eventually 
to be taken. 

We have evaluated (4)-(7) in the replica symmetric 
(RS) Q and the 1st step replica-symmetry breaking 
(IRSB) approximations. In RS one assumes Pa = p for 
the 'polarization', and qaa = Q and qab = q for a ^ b 
for the diagonal and off-diagonal entries of the Edwards- 
Anderson matrix. These are determined from 

P={{V))., C={^^) , q={{vr).. (8) 
\ y/q dz I z 

Here C = /3(g — q), and (. . )z denotes an average over a 
standard Gaussian z while (. . .) is now a Gibbs average 
corresponding to the RS single-site potential 

U,s{v)^-h^sv-]-Cv'' + -G{v) , (9) 
2 7 

with /iRs = Jqp + ^Jqz. The RS approximation thus 
describes a Gaussian ensemble of single-site potentials 
L/cff (w), with parameters p, q and C which are determined 
self-consistently through (8). 

The system described by (8)-(9) exhibits a transition 
from an ergodic phase with p = q = to a glassy phase 
with g 7^ at some temperature Tg depending on Jq 
and 7; see Fig. 1. If Jq is sufl[iciently large, a transition 
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to a polarized phase with p may also occur. The 
assuniption of RS is not always correct. Replica sym- 
metry breaking (RSB) occurs at low temperatures and 
large 7. The location of the instability against RSB is 
given by the AT criterion (3^ { { (v^) - {vf f ) ^. 

In the IRSB approximation, which is expected to consti- 
tute a major step towards the full solution, Ucs{v) is of 
the same form (9), albeit with h^s replaced by /iirsb — 
JoP + ^/go 20 + Vli - 90 zi, and C by C = (3{q - qi), 
where we use standard notation for the entries of the 
q-matrix. Along with a so-called partitioning parame- 
ter m, they are determined from a more complicated set 
of fixed point equations 0,^. Now zq is a Gaussian, 
whereas the distribution of zi is more complicated. 




DWP 



These are barrier height V, asymmetry A, distance d 
between the wells, and the tunnelling matrix element 
Aq — hioo ex.p{— X) , with A = d{2moV/?i'^)^/^ , loq a char- 
acteristic frequency (of the order of the frequency of har- 
monic oscillations in the two wells forming the DWP) and 
mo the effective mass of the tunneling particle. For larger 
I'^rs/irsbIj the C/eff(w) only exhibit single-well forms. 

In contrast to the main assumption of the STM that 
A and A are uncorrelated with P(A, A) ~ Pa, we find 
A and A to be strongly correlated; both are functions of 
one random variable /irs/irsb in RS/IRSB respectively, 
a feature which persists at all finite levels of RSB. By the 
same token, SPM expansion-coefficients of Ues{v) are 
correlated. The correlations between A and A and like- 
wise those between the SPM coefficients can be weakened 
(but not eliminated) by introducing local randomness, 
i.e., by making either 7 or other parameters of the on- 
site potentials z-dependent in a random fashion. It has 
been demonstrated in | [T^ ] that our model remains solv- 
able with these modifications. The RS distributions -P(A) 
and P(A) have been presented in Both have singu- 
larities at their upper boundary, the former an integrable 
divergence, the latter a cusp singularity. 
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FIG. 1. (a): Glass temperature Tg as a function of 7 for 
Jo — 0. (b): T — phase diagram in RS. E denotes the 
T = 0-limit of the ergodic phase, G the glassy phase, and P 
a phase with macroscopic polarization. The full line is the 
AT line. Below the dot-dashed line is the region with DWPs. 
The short-dashed line separates the glassy phase G from the 
phase P with macroscopic polarization. 

Interestingly, the system can also exhibit a collection of 
DWPs in its T = potential energy surface. That is, for 
suitable external parameters, some of the L'efi(i') have 
DWP form. From (2), (9) it is clear that the condition 
for this to occur is 7C > 1; this region is marked DWP 
in Fig. lb. DWPs exist for not too large |/irs/irsbIj 
and the distribution of their characteristic parameters 
derives from the Gaussian nature of z (the joint {zo, zi)- 
distribution in the IRSB case) and can be computed. 



FIG. 2. Low T specific heat in the G-phase for 7 = 10, 
and Jo = in the RS (full line) and IRSB (broken line) 
approximation. Note that the location of the bump at 
Tj, ~ 6 X W~^Tg is reasonable for many glasses. E.g. for 
a-Si02 this would correspond to Tt — 8.6 K. 

Even though P(A, A) is different from what is assumed 
in the STM, we find that the contribution of the tunnel- 
ing excitations to the specific heat exhibits an extended 
range of temperatures where it scales linearly with T. By 
solving the Schrodinger equation of a particle moving in 
the potential Ucs(v) and by computing the specific heat 
beyond the two-level approximation, we obtain a bump 
in a C(T)/T^-plot (Fig. 2), as it is typically also ob- 
served in experiments. Its microscopic origin is related 
to "harmonic" excitations in DWP and SWP structures 
whose frequency depends only weakly on /irs/irsb, thus 
A, and which gives rise to a peak in the density of states 
(DOS). For DWPs with larger asymmetries, the role of 
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tunneling is, indeed, mainly to prevent level crosssings of 
these "harmonic" excitations as A is varied. 

Within our model there is a clear distinction between 
the local G{vi) contributions to J7pot and the global ones 
mediated by the Jij . It is the influence of the latter which 
makes the ensemble of local potential energy configura- 
tions a collective affair and which can thus be regarded 
as a source of universality of glassy low-temperature 
physics. The effect of the on the Ucfi{vi) is to produce 
a linear term — /irs/irsb which creates a broad range of 
asymmetries and thereby a nearly constant contribution 
to the DOS at low energies, and a renormalization of the 
harmonic contribution through —^Cv^ , which is respon- 
sible for the very occurrence of DWPs. Conversely, local 
properties embodied in G{v) are mainly responsible for 
the "harmonic" excitations creating peaks in the DOS on 
top of the unstructured part. These cannot be expected 
to be universal, which fits rather nicely with the experi- 
mental observation that both the bump in the C{T)/T^ 
plot and the corresponding plateau in k{T) (presumably 
created by resonant scattering of phonons from these har- 
monic excitations) do show a stronger material depen- 
dence. It is not unreasonable to assume that such a dis- 
tinction between local and global contributions to Upot 
can be made for real materials as well (perhaps with ran- 
domness in the local contribution), the long-range global 
one being, e.g., of elastic origin. Interestingly, on intro- 
ducing local randomness such as non-degenerate distri- 
butions for the 7, we have obtained slightly super-linear 
low-T specific heats, with exponents ranging between 1 
and approximatly 1.5 [ p^ , which fits nicely with exper- 
imental observations. It is currently under systematic 
study. As mentioned above, local randomness will also 
weaken the perfect correlation between A and A which is 
also desirable in view of the experimental situation Q. 

From the nearly constant DOS at low energies we can 
determine the constant Pq [0- For parameters as in 
Fig. 2, we find Pq - 6- 10-3(fcB A')" Vatom in RS (aprox- 
imately half this value in IRSB), which is roughly three 
orders of magnitude too large for most glasses. While 
this may be regarded as a deficiency, we believe that it 
does not invalidate our approach (no tuning of parame- 
ters whatsoever has yet been attempted) but rather that 
it may be used to guide further detailed modelling: again, 
local randomness turns out to be a possible cure. 

We have proposed and solved a simple model which ex- 
hibits both, an amorphous low-temperature phase, and 
a collection of DWPs and SWPs in its potential energy 
landscape, and which reproduces characteristic glassy 
low temperature anomalies, though, as a model based 
on localized degrees of freedom it cannot be expected to 
present a faithful picture in all details of the ergodic (liq- 
uid) phase. Also, we do not presently regard our model 
as providing a realistic description of specific substances 
but rather as representative for a new and promising way 
of looking at glassy physics — at low temperatures and 



beyond. For instance, models of the type proposed will 
also exhibit interesting dynamical properties at or near 
their glass temperature Tg and possibly throughout their 
glassy phase, which might be investigated, e.g., along the 
lines of |20|. Finally, we believe that our approach does 
hold the potential for eventually including realism; the 
search for good models has only just begun. 

We would like to thank C. Enss, H. Horner, S. Hun- 
klinger, J. van Mourik, P. Neu, O. Terzidis, J. Urmann, 
and A. Wiirger for numerous illuminating discussions. 
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